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We find the Weyl law followed by the eigenvalues of contractive maps. An important property 
is that it is mainly insensitive to the dimension of the corresponding invariant classical set, the 
strange attractor. The usual explanation for the fractal Weyl law emergence in scattering systems 
(i.e., having a projective opening) is based on classical phase space distributions evolved up to the 
quantum to classical correspondence (Ehrenfest) time. In the contractive case this reasoning fails 
to describe it. Instead, we conjecture that the support for this behavior is essentially given by the 
strong non-orthogonality of the eigenvectors of the contractive superoperator. 
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I. INTRODUCTION 

The study of open quantum systems has recently be- 
come a very active field [1]. The reasons are many, includ- 
ing the development of quantum information and com- 
putation 0, Q, quantum optics and scattering systems 
[4, 5]. Particularly in this latter case the fractal Weyl 
law has been proposed. This law predicts the way in 
which the long-lived resonances of these systems grow 
as a function of h. The fundamental ingredient is the 
classical invariant set, which in this kind of systems is 
the repeller, i.e., the set of trajectories non-escaping in 
the past and in the future. In fact, this law says that the 
number of long-lived quasibound states is proportional to 
h~( 1+dH \ where du is the partial Hausdorff dimension of 
the repeller [fj. 

There is a vast literature that has contributed to gain 
confidence on this conjecture by means of numerical tests 
conducted on many systems [7]. However, open quantum 
maps have been the main tool in these studies, as they 
offer great simplicity in the calculations without losing 
much generality |8l-[lQ|. For them, the fractal Weyl law 
predicts that the resonances grow as h~ d , where d is the 
partial fractal dimension of the repeller. But if the way 
to open the system is nonprojective the available litera- 
ture is very scarce. Recently [ll| this situation has been 
analyzed for dissipative quantum operations that can be 
thought as a phase space contraction leading to dissipa- 
tive dynamics [l2[ • In that work a dissipative baker map 
has been studied, where all classical initial conditions 
asymptotically fall on a strange attractor. The quantum 
counterpart has been implemented by means of a noise 
superoperator written in terms of Kraus operators [l3| . 
The number of long-lived resonances has been found to 
behave in a rather different way compared to the usual 
prediction of the fractal Weyl law. In fact, this number 
grows as a power law in ft, but the exponent is mainly 
insensitive to the dimension of the fractal invariant set. 

In this work we analyze this behavior in depth. We 
find the Weyl law for the spectra of contractive noise. In 
order to explain its emergence and discrepancies with the 
usual fractal Weyl law, we first follow the same steps than 
in the case of scattering systems (i.e., having a projec- 



tive opening). This is done in terms of an initial classical 
distribution (that in this case shrinks following the asso- 
ciated dissipation) evolved up to the quantum to classical 
correspondence time Tghr, the Ehrenfest time. We pro- 
pose a theoretical expression for this time based on dy- 
namical considerations and confirm its validity by means 
of the exploration of the classical phase space distribu- 
tions and the eigenvectors of the contractive superoper- 
ator. However, this reasoning does not lead to a satis- 
factory explanation. We conjecture that the strong non- 
orthogonality of the right eigenvectors is the main reason 
behind this behavior. 

This paper is organized as follows: in Section [II] we 
briefly describe the dissipative model that we have used 
and give the expression of the Weyl law for the contrac- 
tive baker map. In Section Hill the numerical results are 
analyzed and we explore possible explanations for the 
emergence of the here obtained Weyl law supported by 
studies of the phase space distributions and the proper- 
ties of eigenvectors. Finally, we give our conclusions in 
Section \N\ 



II. THE WEYL LAW FOR CONTRACTIVE 
MAPS 

As in our previous work [ll| we have investigated the 
spectral behavior of the dissipative baker map, which is 
defined on the 2-torus T 2 = [0, 1) x [0, 1) by 



B(q,p) 



(2q,ep/2) 
(2q-l,(ep- 



if < q < 1/2 
l)/2) if 1/2 < q < 1 



(1) 



Besides contracting the torus in the p direction by a 
e factor, this map stretches the unit square by a factor 
of two in the q direction, squeezes it by the same factor 
in the p direction, and then stacks the right half onto 
the left one. As a result a strange attractor sets in after 
a few time steps regardless of the nature of the initial 
condition. 

The first step to quantize it is to impose on any 
state periodic boundary conditions on the torus, for 
both the position and momentum representations. Then, 
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we take (q + = e* 27rx « (q\i/j), and (p+l|^) = 

e l2irx p (p\ip), with Xg, Xp £ [0, 1). There is a finite dimen- 
sion N = (2t:K)~ 1 for the corresponding Hilbert space 
and a discrete set of position and momentum eigenstates, 
which is given by \q 3 ) = \{j + Xq)/N) (j = 0, 1, . . . iV-1), 
and \p k ) = \(k + Xp)/N) (k = 0, 1, . . . TV - 1), whose 
eigenvalues are q^ p^. A discrete Fourier transform, i.e. 



(G N )kj = (Pk\qj) 



1 / -i2ir, W7 s 



relates these sets. We take anti-symmetric boundary con- 
ditions, this meaning Xq = Xp = 1/2- For an even TV- 
dimensional Hilbert space, the quantum baker map is 
defined in the momentum representation as [3, [H| 
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with Bn a unitary matrix (closed quantum baker map). 

We introduce dissipation by means of a non-unital 
quantum operation [12| implemented by an iV 2 x iV 2 
Kraus superoperator of the form: 



JV-l 



(3) 



FIG. 1: (Color online) Weyl law for contractive noise: 
/long-lived as a function of (y C ute)/N. Results for e = 
0.8, 0.7, 0.6 and 0.4 are represented by means of up triangles 
(in red), down triangles (in black), dots (in blue) and squares 
(in magenta), respectively. 



Quantum operations act on the density matrix, de- 
notes the place where this later must be inserted in or- 
der to implement the corresponding quantum operation. 
Here 



N-l 



-M(l- e )M| pi _ M )( pi | ( 4 ) 



are operators that induce transitions towards the momen- 
tum state \pi = o).The coupling constant e has the same 
value as the dissipation parameter of the corresponding 
classical map. M is a trace preserving (52 ^ Aj^A^ = 1) 
and non-unital (52 ^ A^A^ ^ 1) superoperator, which de- 
scribes a process contracting phase space volume. The 
complete quantum dissipative dynamics is obtained by 
composing M with the unitary map ([2]), 



$ = (B N ® B* N ) o M. 



(5) 



In this work we have computed the eigenvalue spec- 
trum of superoperoperator (|5]) for different values of the 
contraction parameter (e = 0.8, 0.7, 0.6, 0.4) and of 
the dimension (90 < N < 180). For each case we 
have counted the number of complex eigenvalues A (with 
|A| = exp(^)) with a decay rate 7 smaller than a given 
value 7cut- The data are collected in Fig. [1] which dis- 
plays the fraction of resonances fiong- lived as a func- 
tion of e, N and the cut-off value 7 cut (in a wide range 
2<7cut<14). 

By fitting these numerical results, we obtain a remark- 
ably compact and simple expression: 



/long— lived (e,7cut, N) 



N, 



7<7o 



N 2 



C (e 7 cut) 2i/ (iV 2 )- 



(6) 

The values of C and v, for four different values of e, 
are given in Table HI In the fourth column we display 
the semiclassical prediction ^ sc , which will be analyzed 
in Section HTT1 



TABLE I: Values of the fitted coefficients C (column two) 
and v (column three) for different values of e. The fourth 
column displays the semiclassical prediction is sc described in 
Section [nil 



e = 0.8 


C = 5.3 


v = 0.72 


v S c = 0.24 


e = 0.7 


C = 4.4 


v = 0.76 


vsc = 0.34 


e = 0.6 


C = 4.7 


v = 0.79 


vsc = 0.42 


e = 0.4 


C = 5.3 


v = 0.85 


v S c — 0.57 



These findings generalize the ones obtained in On 
the one hand they confirm the existence of a power law 
dependence of Ny ^p ut on N with an exponent which, in 
a meaningful range of validity, is fairly insensitive to the 
value of the dissipation parameter e. On the other hand, 
they hint (within a precision of 20%) on a very simple 
dependence of the prefactor with both e and the cut-off 
value 7 cut . We will leave the analysis of this prefactor, 
which is in general believed to be system-dependent, for 
future work [18| and concentrate in the following on the 
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scaling of 7 ^2 cut with N. We will seek for an expression 
of in order to determine to which extent this expo- 
nent can be related to the underlying classical dynamics. 
For this we will follow an approach analogous to the one 
used in the formulation of the fractal Weyl law for chaotic 
maps with a projective opening [loj and discuss its limi- 
tations in the case of a contractive noise. 



III. CLASSICAL AND QUANTUM SUPPORT 
FOR THE EIGENVALUE STATISTICS 

A heuristic formulation of the fractal Weyl law for 
chaotic maps with a projective opening is based on the as- 
sumption that the number of long-lived resonances (asso- 
ciated with the classical repeller) scales as the volume of 
the evolved initial classical distribution up to the Ehren- 
fest time, that is, the volume of a finite (Ehrenfest)time 
repeller [10[. This volume can be calculated by a combi- 
nation of two exponential laws that relate the probability 
to reside in the system (non escaping trajectories) and 
the quantum to classical correspondence. 

In the case of a contractive noise the connection be- 
tween the long-lived resonances and the structure of the 
classical invariant also exists. In particular we have veri- 
fied in that the Husimi representation of the projec- 
tor corresponding to the eigenfunctions with slow escape 
rate concentrates on the phase space region correspond- 
ing to the classical strange attractor. It seems then nat- 
ural to generalize the considerations usually applied to 
chaotic maps with a projective opening to the contrac- 
tive case and investigate whether this scheme succeeds in 
accounting for the Weyl law o f eq .fl6]). Our starting point 
will be the following relation jl6j : 

/long-lived (e, 7cut, N) ~ Aj? las (7) 

between the fraction of long-lived resonances and the vol- 
ume of the attractor A c \ as which shrinks exponentially 
until the Ehrenfest time according to: 

^cias = exp -(7d T Eh r). (8) 

Notice in eq.© the square (instead of linear) depen- 
dence on A c \ as , which is due to the use of the superoper- 
operator formalism to model the contractive noise. The 
classical decay rate j c \ and the correspondence (Ehren- 
fest) time TEhr are then the two main ingredients of this 
approach that should be evaluated. 

The classical decay rate can be easily calculated by 
following the time evolution of a uniform distribution in 
phase space under the action of dissipation. It is straight- 
forward to see that after t time steps the original dis- 
tribution will occupy 2 t fringes in the q direction, each 
fringe having a width (-)*. Hence the total phase space 
area occupied by the distribution as a function of time is 
^-cias = e _t7cl = e*, and then the classical decay rate is 
given by j c \ = -Inc. 



Determination of the Ehrenfest time is a more subtle 
issue. Understood as the time at which the quantum and 
the classical descriptions differ, we can start our reason- 
ing following the lines of what is done in the case of area 
preserving maps. In fact, there are two different ways to 
conceive this correspondence time. The first one is the 
time Tehri at which a given initial semiclassical distribu- 
tion (a coherent state of width y/h, for instance) spreads 
up to the border of the system along the unstable direc- 
tion (manifold). This time is related to the expansive 
Lyapunov exponent Ai, such that Tshri °c ^x^- On the 
other hand, the time TEhr2 is that corresponding to the 
initial distribution shrinking along the stable direction to 
a size of the order of the Planck cell (1/N). This time 
is related to the contractive Lyapunov exponent A2, such 
that TEhr2 oc j^- Of course, in the case of an area 
preserving map Ai + A2 =0 and both times coincide. 
However, under a contractive noise, our dissipative map 
gives Ai = ln2 while A2 = — . Hence, we propose 
the shortest TEhr2 oc InN/ln^ as the global quantum to 
classical correspondence time for this map. 

In order to verify this assumption, we have numerically 
estimated the correspondence time. This can be accom- 
plished quite easily by evaluating the overlap O c i-q be- 
tween the finite time classical attractor and the Husimi 
distribution of a uniform initial state evolved up to the 
same time. If we exploit the fact that for the baker map 
the interesting features of the distribution (namely its 
fractality) are only in the p coordinate we can notably 
simplify this task. In fact, we just calculate the norm of 
the evolved wavefunction, restricted to the region occu- 
pied by the classical distribution at any given time. As a 
result we have obtained Fig. [2] where these overlaps are 
shown as a function of the map iterations. We have found 
that, besides small fluctuations and the lack of precision 
inherent to the discrete time steps of the map, the results 
confirm our theoretical prediction (see vertical lines as a 
guide). 

Inserting the expressions of j c \ and TEhr2 in eq.flSJ) 
gives A c i as = N~ Usc , where z/ sc = 2 — d, and d = 
1 + In (2) /(In (2) — In (e)) is the fractal dimension of the 
classical attractor. The values of the semiclassical v sc 
are listed in the fourth column of Table HI showing a dra- 
matic discrepancy with the values obtained by fitting our 
numerical results with eq. ([6]). Besides an overall factor 
of ~ 2 between v and ^ sc , the semiclassical exponent 
shows a dependence on e (via the fractal dimension of 
the attractor) which is absent in the fitted v which are 
practically constant. Then, it becomes clear that the way 
of reasoning that has provided with a reasonable expla- 
nation for the emergence of the usual fractal Weyl law 
for systems subjected to projective noise can no longer 
be applied to contractive dynamics. We are now faced 
with the question of where this discrepancy comes from. 

At the basis of eq. ([7j) is the assumption that the num- 
ber of long-lived quantum states can be approximated 
by the number of Planck cells which fit into the phase 
space volume of the classical invariant set. This, in turn, 
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FIG. 2: (color online) Overlap between the phase space region 
occupied by the quantum and classical attractors as a function 
of time t (map iterations). Upper panel corresponds to e = 
0.4, lower panel to e = 0.8. Results for N = 100, 200, and 400 
are represented with squares (in black), dots (in blue) and up 
triangles (in magenta), respectively. Vertical lines show the 
corresponding theoretical values of Tem for N = 100, 200, 
and 400. (full, dashed and dotted lines, respectively). 



supposes that to a good approximation the eigenfunc- 
tions supported by this set are non-overlapping. Even 
though we cannot strictly speak of orthogonality, since 
the operators describing open systems are not normal, 
we know that in the case of projective openings the 
long-lived eigenfunctions are quasi-orthogonal (while the 
short-lived ones present a high degree of degeneracy). 
This explains the success of the fractal Weyl law in the 
projective case. In the case of contractive dynamics we 
will investigate this point by defining the overlap ma- 
trix Pij = Tr(R\Rj), where Ri are the right eigenstates 
corresponding to the superoperator $ of Eq. ([5]) (this is 
not to be confused with the biorthogonality of the right 
and left eigenfunctions of a superoperator, which states 
that Tr(L[Rj) = Sij). The overlap matrix elements cor- 
responding to the contractive map with N = 180 and 
e = 0.4, 0.6, 0.8 for the 200 longest-lived eigenstates 



are displayed in panels (a), (b), and (c) of Fig. [3j respec- 
tively. A grayscale is used to represent them , going from 
white corresponding to value to black corresponding to 
the maximum values. We observe that the off-diagonal 
elements are clearly non negligible. Moreover, their value 
grows with the contractive power of the corresponding 
map (as e decreases). For comparison we show in panel 
(d) the overlap matrix for a projective case, obtained by 
opening the baker map along two symmetric bands in the 
q-direction, of width Sp = 0.1 and centered at p = and 
p = N — 1. In this case, as expected, the matrix is almost 
diagonal. 




FIG. 3: Overlap matrices Pij of the first 200 right eigenstates 
with N = 180 (ordered by decreasing modulus of the eigen- 
values) corresponding to the contractive map for e = 0.4, 0.6, 
and 0.8, (panels (a), (b), and (c) Respectively). For compar- 
ison we show the same overlap matrix but for a projective 
opening that amounts to 0.2 of the phase space. Only the 
upper half of the matrices is shown. 



The different degree of non-orthogonality of the long- 
lived resonances in both models is also reflected in the 
phase space distribution of these states. In panel a) and 
c) of Fig. [H we show the sum up to 7 cut of the Husimi 
representation of the longest-lived right eigenstates: 



g {z\R,\z)(z\R\\z) 

7=0 



(R-y | R\ ) 



(9) 



with {z\R-y\z) = Tr(R\,\z){z\) where \z) are coherent 
states centered at z = (q,p). Panels b) and d) display 
the analogous sum (|9|) corresponding to the Husimi repre- 
sentation but of the Schur eigenvectors, which constitute 
the orthogonal basis associated with the eigenvalues A 
with |A| > exp(^jf^). 

In the case of the contractive map (upper line) we ob- 
serve that the area of phase space occupied by the sum 
of the Husimi distributions is smaller than the area cor- 
responding to the subspace spanned by the Schur decom- 
position. This is a clear sign of the non-orthogonality of 
the eigenstates for this kind of superoperators [17j. 
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On the contrary, the lower panels (c) and (d) show that 
for the case of a projective opening both distributions 
look much the same, indicating that the assumption of 
quasi-orthogonality for the long-lived eigenfunctions is 
justified. 



(a) 



(b) 
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FIG. 4: In panel (a) we show the sum J9J of the Husimi repre- 
sentation corresponding to the first 600 right-right eigenvec- 
tors for the contractive map at e = 0.6 for N — 180. Panel 
(b) displays the analogous sum corresponding to the Husimi 
representation of the Schur eigenvectors. For comparison, in 
the lower panels (c) and (d) we show the same distributions 
than in panels (a) and (b) respectively but for the projective 
opening case of Fig. [3] (d). 



pie dependence on (7 cu te)/iV has also been obtained for 
a dissipative kicked top map on the sphere. We were not 
able to explain the emergence of this law by means of 
the usual line of reasoning applied to the projective case. 
Very simply put, the idea is counting resonances. This 
has been traditionally accomplished by partitioning the 
phase space volume occupied by a finite time classical 
invariant set (the repeller). In fact, implies a pseudo or- 
thogonality of the long-lived eigenstates. We could verify 
that this is indeed the case for the projectively opened 
baker map, a system that has been paradigmatically used 
in the fractal Weyl law literature. But when it comes 
to the dissipative baker map used in this work, we have 
clearly identified a high degree of non-orthogonality. This 
is the main reason behind the failure of the usual reason- 
ing for explaining the emergence of the Weyl law. 

As a result, we think that a new method for counting 
the long-lived resonance other than just partitioning the 
corresponding volume in phase space into Planck cells, is 
the key to understand the statistical behavior of contrac- 
tive maps. In the future, we hope to find a theoretical 
explanation for it, including the one of the prefactor and 
the dependence on e and 7 cut Isj. 
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